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Abstract 

We study a non local version of the sine-Gordon model connected to a many- 
body system with backward and umklapp scattering processes. Using renormal- 
ization group methods we derive the flow equations for the couplings and show 
how non locality affects the gap in the spectrum of charge-density excitations. We 
compare our results with previous predictions obtained through the self-consistent 
harmonic approximation. 
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1 Introduction 



In general, Quantum Field Theories have been built in the context of local models. 
However, there exist physical situations that lead to non local interactions in a straight- 
forward way. Let us mention, for instance, Wheeler and Feynman's description of charged 
particles [JJ, string theories with non local vertices |2J, and non local kinetic terms that 
appear when bosonizing fermions in (2+1) dimensions jl]- As we shall see, several 
recently considered non local field theories are related to the study of electronic systems 
in one spatial dimension (ID) jH] Indeed, in recent years the physics of ID systems 
of strongly correlated particles has become a very interesting subject since one can take 
advantage of the simplicity of the models at hand and, at the same time, expect to make 
contact with experiments. For instance, the recently discovered carbon nanotubes are 
perfect experimental realizations of ID conductors [Zj. On the other hand, as the dimen- 
sionality of a system decreases, charge screening effects become less important and the 
long-range interaction between electrons is expected to play a central role in determining 
the properties of the system. In fact, from a theoretical point of view the effects of long- 
range interactions have been recently discussed in connection to several problems such 
as the Fermi-edge singularity jH] , the insulator-metal transition [0] , the role of the lattice 
through umklapp scattering and size dependent effects [10], etc. In the specific context 
of carbon nanotubes, a low-energy theory including Coulomb interactions has been also 
recently derived and analyzed ^JJ. Non local fermionic models have been also used in 
the study of fluctuation effects in low-dimensional Spin-Peierls systems [T2"] . 

As shown in |13j . starting from a non-local and non-covariant version of the Thirring 
model P3], in which the fermionic densities and currents are coupled through bilocal, 
distance-dependent potentials, one can make direct contact with the "g-ology" model 
[To] currently used to describe different scattering processes characterized by coupling 
functions g\, g2, g% and g&. When one bosonizes this theory by either operational or 
functional methods, due to the contributions of g\ (backscattering) [TE] and g 3 (umklapp) 
[T7j one finds an even more drastic departure from the local case. Indeed, instead of the 
well-known integrable sine-Gordon model (SG) one gets a non local extension of it, which, 
as far as we know, is not exactly solvable. Recently, in ref.[T3|, the physical content of 
this model was explored by using the self-consistent harmonic approximation (SCHA) 
|18j . As it is well-known the SCHA is a non controlled approximation, i.e. there is 
no perturbative parameter involved. It is then desirable to have an alternative analysis 
of this problem. This is the main motivation of the present work. We will apply the 
renormalization group (RG) technique ^H[, usually employed in local cases, to the non 
local sine-Gordon model (NLSG) mentioned above. For simplicity we shall assume that 
non locality plays a role only in umklapp interactions (g^) whereas all other potentials 
are local, i.e. proportional to delta functions. In Section 2 we briefly show how the 
NLSG action is obtained from the non local Thirring model. In Section 3 we derive the 
RG equations and compute the gap in the charge- density spectrum. This allows us to 
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determine the effect of non contact g 3 couplings. Finally, in Section 4, we discuss our 
results. 



2 The Model 

Let us sketch the derivation of the NLSG action. We start from the fermionic (1+1)- 
dimensional Quantum Field Theory with Euclidean action given by 

S = So + Sf s + Sbs + S us (2-1) 

where 

S = J d 2 x^i$^ (2.2) 

is the unperturbed action associated to a linearized free dispersion relation. The contri- 
butions of the different scattering processes can be written as 

Sfs = -y J d 2 xd 2 y(^)(x) V w (x,y) (*7„*)(v) (2-3) 

and 

S bs + S us = ~YJ Sxd2 y r ^)( x ) Uv)(x,y)(* T^)(y) (2.4) 

where the j^s are the usual two-dimensional Dirac matrices and r = 1, Ti = 75. The 
coupling potentials V( M ) and are assumed to depend on the distance | x — y | and can 
be expressed in terms of "g-ology" parameters as 





(x,y) = 


-^(92 + 94}(x,y) 


V(i) 


(x,y) = 


-2(92 -g*)(x,y) 


U(o) 


(x,y) = 


■^2(93 + 9i){x,y) 
■^2(93 -9i)(x,y) 


U(i) 


(x,y) = 



(2.5) 
(2.6) 
(2.7) 
(2.8) 

In the above equations g and g' are just numerical constants that could be set equal 
to one. We keep them to facilitate comparison of our results with those corresponding 
to the usual Thirring model. Indeed, this case is obtained by choosing g' — and 
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V(o)(x,y) = Vm(x,y) = 5 (x — y). On the other hand, the non-covariant limit g' = 0, 
Vm(x,y) = gives one version ((g 2 = g^) of the TL model [2Uj . 

The terms in the action containing g% and g^ represent forward scattering events, in 
which the associated momentum transfer is small. In the gi processes the two branches 
(left and right-moving particles) are coupled, whereas in the g± processes all four partic- 
ipating electrons belong to the same branch. On the other hand, g% and g^ are related 
to scattering diagrams with larger momentum transfers of the order of 2k p (bs) and 4k p 
(us) respectively (this last contribution is important only if the band is half-filled). For 
simplicity, throughout this paper we will consider spinless electrons. The extension of 
our results to the spin-1/2 case with spin-flipping interactions, though not trivial, could 
be done by following the lines of ref. fU\ . 

At this point we consider the partition function Z expressed as a functional integral 
over fermionic variables. The implementation of a generalized Hubbard-Stratonovich 
identity 6\ allows to write Z in terms of a fermionic determinant. Although this determi- 
nant is highly non trivial, one can combine a chiral change in the fermionic path-integral 
measure with a formal expansion in g' in order to obtain a bosonic representation (See 
ref. [T3| for details). One thus establishes an equivalence between the original fermionic 
action and the following bosonic action depending on five scalars $, rj, C , C\ and ip: 



2 

/ p3i [*(p)*(-!>M(!>) + l(pH-p)B(p) + HpH-p)C(p) + <fi(p)tfi(-p)j} 



+ 2 
where 



and 



d 2 x d 2 y C^x)!!, J(x,j/)C M (y) + 



A(p) 
B{p) 



g'Ac 



71 



d 2 x [C (0) (x)fo(x) + iC {l) {x)h{x)} (2.9) 



plv m 1 (p)+pi(v (0) 1 (p) + 9 -) 
pIv {0) 1 (p)+pI(v w 1 (p)- 9 -) 



C(p)=p p 1 lV {0] 1 (p)-V [ - ) 1 (p) + y - 



fo(x) 
fi(x) 



cos((v47T(/9 - 2ig$)(x)) 
sin((A/47ry? - 2ig<f>)(x)) 



(2.10) 
(2.11) 
(2.12) 

(2.13) 
(2.14) 



Since the integrals in C and G\ are quadratic these fields are easily integrated out 
and one gets 



Z = M j V<S>Vr]Vipe 



-s, 



eff 



(2.15) 
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with 

S eff [$,ri,<p]=S + S int (2.16) 

where 



/d 2 2 
^2 [^(P)^(-P)^(P)+»7(P)»7(-P)5(P)+^(P)»7(-P)C(P)+¥'(P)¥'(-P)^-] (2-17) 



(Ac) 



d 2 xd 2 yg x {x,y) cos ^V47r(y?(:r) - <p(y)) -2ig($(x) - $(y)) 
d 2 xd 2 yg 3 (x,y) cos \V^((p(x) + <p{y)) - 2ig(${x) + $(y)) . (2.18) 



2tt 2 

(Ac) 2 
2tt 2 

It is now convenient to diagonalize the quadratic part of the effective action by intro- 
ducing the fields (, x an d 4>'- 



2igBp 2 



g A + 2Bg 2 p 2 







7/ 



-iC 



C 



igCp 2 



1 



2£(? A + 2Bg 2 p 2T g 

-c + A 



where we have defined g 2 = g 2 / 71 " an d = CO ) 2 — 4A(p)_B(p). We then obtain 



(2.19) 
(2.20) 
(2.21) 



. _ 1 f d 2 i 
'°-2j (27 



d 2 p_ 

? 



C(p) p 2 + 



25^ 



/,2 



C(-p) + x(p)^x(-p) + 0Cp) A + 2£? ~ 2p 



(2.22) 



2tt 2 



y d 2 xd 2 yg 1 (x,y) cos v 7 !^^) - 0(y)] 
(Ac) 2 



2vr 2 



d 2 xd 2 y^ 3 (^y)cosV47r[0(x) + 0(y)]. (2.23) 



One can see that the ( and x fields become completely decoupled from <fi. Moreover, it 
becomes apparent that the ^-dependent piece of the action Si n t is the only one containing 
potentially relevant contributions (i.e. gapped modes). 
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3 RG treatment of the non-local umklapp coupling 



In this Section we shall focus our attention on the non local action derived above. For 
simplicity, from now on we will consider the case in which (72(f) and (74(f) are constants 
(local forward scattering) and g\ = 0, g 3 7^ 0, i.e. a pure non local umklapp interaction. 
Thus, we start from the action 

S[4>] = J 7^*(P)^#-P)-^ / d 2 xd 2 yg 3 (x,y)cosV^{^x) + <P(y)] (3.1) 
with 




(3.2) 
(3.3) 



l + ^ + ^fl + ---)> (3-4) 



7T 7T / V 7T 7T 



where we have now expressed all formulae in terms of g coupling functions. In the local 
case ((73 = 5 2 (x — y)) the action (j3.1j) corresponds to the well-known sine-Gordon model, 
which is an integrable, exactly solvable field theory. In particular, a RG analysis shows 
that the "stiffness constant", K has to be lower than .5 in order to have a relevant cosine 
interaction, i.e. to have a gap in the spectrum. Recently, by reinterpreting Bethe-Ansatz 
results, Zamolodchikov obtained the exact expression for this gap [221 • Unfortunately, as 
far as we know, the present non local version of the theory is not exactly soluble and one 
is then forced to consider an approximation. In ref. ^Hj a SCHA was employed in order 
to estimate the energy gap. But this is a non controlled, non perturbative approximation. 
Besides this general disadvantage, the implementation of the SCHA technique led to a set 
of coupled algebraic equations that could be numerically solved only for very weak non 
locality. It is then natural to try another method to attack the problem and eventually 
improve the approximation. Let us consider the Wilsonian approach to the RG (See for 
instance [El)- First of all we restrict our analysis to a non local interaction of the form: 

9z{x-V)=9a{x\-Vx)K x Q-Vo)> ( 3 - 5 ) 

with 

g 3 (xx) = A <5(xi) - -j^<%5(xi) (3.6) 

where df5 2 (x) is the second derivative of the delta function with respect to x\. At 
this point it is worth mentioning that RG calculations involving fermionic non local 
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interactions already exist in the literature. For instance, the authors of ref. [12] used a 
Jordan- Wigner fermion representation for the ID Heisenberg-Ising model which includes 
not only a non local fermion-fermion interaction but also a linear fermion-lattice coupling. 
Due to the curvature of the band, the non trivial fermion-lattice coupling and the presence 
of both forward and umklapp scattering, comparison of their RG equations with ours is 
not straightforward. Let us point out, however, that in jT2] fermions corresponding to 
nearest-neighbor lattice sites interact through a potential of the form g(q) — C\ X cosg, 
i.e. a unique constant C\ is associated to g(q)- As shown in the above equation ()3.6|) . in 
the present work we are interested in a coupling which depends at least on two constants 
A and e. In fact, the derivation of RG equations for these coupling constants will be our 
next task. 

In condensed matter problems one is usually interested in the physics at long distances, 
compared to a lattice spacing of the order of A -1 . Since, in momentum space, this 
corresponds to small k\ = k, it is natural to consider correlations between fields with 
momenta < k < A/s, with s very large. These are the so called "slow modes" 0<. On 
the other hand, the "fast modes" 0> are those carrying momenta that satisfy A/s < k < 
A. In the present approach to RG these fast modes are integrated in the path integral 
framework, giving rise to an effective theory depending only on slow modes. From this 
action, in turn, one can read the flow equations for the couplings. Indeed, writing the 
initial action as So — Si n t, to first order, after a suitable rescaling of coordinates and 
momenta and a redefinition of the fields (see the Appendix for details) we obtain the 
following relationship between the original and RG transformed actions: 

S int = ^£ J d 2 xd 2 y f[X + A (2 - AK)t - 2Ke t}6 2 (x - y) - - - d 2 x 8 2 {x - y)j x 

x cos V / 47r(0(x) + 4>{y)) = 

(Ac) 2 

2tt 2 J a V"" v ~ A 2 

As usual, imposing the invariance of the action under RG, we get the flow equations 
for the couplings A and e: 

^ = (2 - 4K)X - 2Ke (3.8) 

Tt = < 3 - 9 > 

where Ins = t, with the initial conditions A(0) = Ao and e(0) = eo- The solution of this 
system is elementary, yielding: 

A(t) = (A - Ke )exp[(2 - AK)t] + Ke exp [-AKt] , (3.10) 



d 2 x d 2 y XS 2 (x -y)-^ df5 2 (x - y) cos V4tt(0(x) + <f>{y)) = S int . (3.7) 
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e(t) = e exp[-4JCt]. (3.11) 

From the last equation one clearly sees that the non local piece of the interaction is 
irrelevant, as expected. Concerning the local interaction one sees that it is relevant for 
K < 1/2, i.e. non locality does not modify this well-known condition already found for 
the local SG model. Therefore, for K < 1/2, A will grow with increasing t and there 
will be a gap in the CD spectrum which can be estimated by determining the value t for 
which A = 1. From now on we shall restrict our analysis to the case K < 1/2. The gap 
is then given by AE = A/ie -i , where /i = c/(7r\/2) ~ .198. It is also convenient to define 
the dimensionless gap m = AE/A. From equation (|3.1()jl one can thus derive an equation 
for m which gives the behavior of the energy gap as function of the forward scattering 
potentials (K) and the non local contribution of the umklapp scattering (eo). Before 
analyzing this non trivial equation it seems reasonable to check if it predicts sensible 
results for the local case. To this end we set e = and A = 1 in ()3.1()j) . obtaining: 



m = n\ 2 - AK . (3.12) 

This result can be compared with the exact solution obtained by Zamolodchikov [22] 
and with the approximated result given by the SCHA method jTHj. The corresponding 
expressions for the gap are respectively given by: 

2™ . 1 r(jjy / T{l-2K) 



m scha = -= (4^)™ ii— A 2 - 4 *. (3.14) 
V K 



and 



In order to compare m and m sc h a with mz in an efficient and easy to visualize way, we 
have computed the relative error Am/mz as function of K for both approximations. The 
results are depicted in Fig.l where one sees that our RG computation gives values of the 
gap closer to the exact values for a wide range of the stiffness constant K. Interestingly, 
the SCHA result works well when one approaches the end points of the interval. 

Going back to the case eo 7^ 0, by combining eqs. ()3.8j) and ()3.9|) one readily gets a 
phase diagram in the e — A plane (See Fig. 2). There is a critical line given by A = K e. If 
the initial parameters are tuned to lie on this line, the system will flow to the Tomonaga- 
Luttinger fixed point, at the origin. In this case, of course, the system remains gapless. On 
the other hand, for initial conditions outside the critical line, the system flows to strong 
coupling, giving rise to a gap m, as mentioned above. For simplicity let us consider the 
case A > e and define the variable x = m/m . The gap equation can then be written as 

* 2 --* 2 - 4 * + -(l-^) = 0, (3.15) 
eo e o Ao 
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where v = Xq~ 2K / K. This equation is one of our main results. For fixed values of K 
and A it gives the behaviour of the energy gap as function of the non local contribution 
to umklapp scattering, associated to non contact interactions. The form of this formula 
suggests that it could be easier to handle the inverted equation: 

1 _ X 2 ' AK 

e Q (x) = e cri t v — — -, (3.16) 

^ ^crit % 

where e cr u = Xq/K. In Fig. 3 we show the numerical solution of this last equation for 
Ao = .5 and K = .25. We see that x decreases for increasing eo, in qualitative agreement 
with the SCHA prediction [T3*] . A quantitative comparison between both approximations 
is given in Fig. 4, for the same fixed values of A and K. Since the SCHA result obtained 
in [T3] is valid for small values of the coefficient associated to non locality, we have plotted 
the solutions in the interval < e < 0.1. We find that the gap decay predicted by RG, 
for increasing non locality, is much slower than the one obtained through the Gaussian 
approximation. 

As a final comment, we note that the RG treatment for non local interactions depicted 
in this work can be extended to a more general coupling function including an arbitrary 
number of even powers of p\. In coordinate space such an interaction can be written as 



9s{x) = D- 1 )^ 2 ^) ( 3 - 17 ) 

i=0 



where Ai = eo- The corresponding set of RG equations for the A^s can be obtained from 
the expansion of the cosine integral function which appears when one integrates the fast 
modes after the mode separation (See Appendix). The result is 



^ = (2-4K)Ao-2KX> 



a -i = (2-2j-±K)\ j -2K A, 



^ = (2-2n-AK)X n . (3.18) 



The general solution of this system can be expressed as a combination of exponentials 
and the computation of the gap cannot be done, in principle, in an analytical way. It 
is then illustrative to consider a particular case in which this calculation is simplified. 
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Indeed, for small £, it can be proved by induction that the solutions of this system of 
equations are of the form: 



9K 

X 3 (t) = A;exp[t(2 - 2j - 4K - -p- J2 \ ? )]- (3-19) 

3 i=j+l 



In this limit one obtains the following expression for the gap: 



i 



2-4if-^K- V n ' A° 
^0 ^t=x 



m = /xA , (3.20) 



which is consistent with the conditions Aq ~ 1 and K <C n- 



4 Conclusions 

In this paper we have considered a non local extension of the sine-Gordon model. This 
theory is obtained when one bosonizes a non local and non covariant version of the 
Thirring model used to describe certain Id many-body systems. Since the integrability 
of this non local sine-Gordon model has not been proved, one needs to implement some 
approximation to study its physical content. We have performed a RG calculation up to 
first order in the coupling function g$, which in a condensed matter context is associated 
to the so called umklapp scattering. We obtained an expression for the energy gap as 
function of the non local piece of the interaction e . For purely local interactions (the 
exactly solvable SG) our result seems to be a sensible approximation, improving the 
SCHA predictions for a wide range of forward interactions. In the non local case, in 
which no exact answer is known, we predict decreasing values for the gap for increasing 
values of eo, in qualitative agreement with a previous SCHA computation. We were able 
to give a precise comparison between both approximations in the interval < eo < 0.1, 
showing that the gap decrease, for increasing non locality, is much weaker according to 
the RG computation. Since, as is well-known, the SCHA method is not a controlled 
approximation, the present results contribute to a better understanding of the physics of 
non local field theories. We think that our results are also of interest in the context of 
ID many-body systems (Luttinger liquids) in which most of the previous investigations 
involving umklapp scattering do not consider non local effects associated to long range 
interactions j2H|. 
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5 Appendix 



In order to illustrate the computation leading to the flow equations we first define the 
free bosonic propagator: 

< mm >»- j-J^ mm e-^ = j- ^ (s .i) 

with 

= y £>0 exp[- | ($(p) F(p) </>(-p) - J(p) 0(p))]- (5.2) 

The result is 

5 2 (V> + a) 47T 2 i^f 
< 0(P) 0(g) >o= ^t? 2 2 ■ ( 5 - 3 ) 

Po + w Pi 

The next step is the analysis of S int , as given by the second term of equation (|3.1|) . 
For simplicity, in this Appendix we disregard the overall constant (Ac) 2 / {2tt 2 ) . Going 
to momentum space and performing the separation in slow and fast modes 0< and </>>, 
according to: 

^)=(2^2 / d 2 p<P( P )e^ (5.4) 
0(P) = 0<(p),|Pi| < A/s (5.5) 
0(p)=0>(p),A/s<| Pl |<A, (5.6) 

we obtain 



Sm = J rf 2 xrfyi5f 3 (xi - ?/i)(cos[^^ y d 2 p^<(p) /(p,x,yi)] 



cos [?^3i y rf2 P0>(p)/(p,a; ! l/i)] -sin[-^^ / d 2 p 0<(p) /(p, 



sin [77T~T7 / rf2 P0>(p)/(p,a;,Z/i)]) (5.7) 



(27T) 2 

where 

/(p, x, yi) = e ip - x + e i(P0X0+PlXl) . (5.8) 

Now we expand the functional integral up to first order in #3 and integrate over the 
fast modes. One finds the following results: 

y/An f 

<sin[^-py ^p^p) f(p,x, yi )] > = (5.9) 
10 



and 



< cos 



[j^2 J d 2 p<t>>(p) f(p,x,yi)} >o = 

( v { a (at \ at \\ 1 + C0S (Pi( x i -Vi)) 
= exp l-K J dpi (0(px) - 0{-pi)) 



= exp — K 



2lns + 2Ci(A(x 1 -y 1 ))-2Ci(-(x 1 -y 1 )) 



(5.10) 



where Ci(x) is the cosine integral function and the free vacuum expectation values are, 
of course, taken with respect to fast modes. Rescaling momenta, coordinates and fields 
such that the free piece of the action So remains invariant: 

p = S p, x = s^x, <f>'(p) = s~ 2 0<(p' '/s), (5.11) 

and using the fact that In s — t, s ~ 1 + 1, S int can be written as 

S in t = J d 2 x' dy[g 3 (s(x[-y[)) cos[^0^ J d 2 p<p < (p) f(p,x,yi)]x 

x (1 + [(3 - AK) + KiAix, - y[)) 2 - ^(A(x[ - y[))*]t), (5.12) 

where we have used the power expansion of the function Ci(x). Finally, using the explicit 
expression for g 3 in terms of Aq and eo we obtain 



S int = j d 2 x d 2 y'[(X + A [2 - 4K]t - 2Keot)5 2 (x[ - y[)- 

-9?<$Vi - y'j] cos y/^[(f>(x) + 4>{y')] 



eo-AKe t 2 , 



A 2 



= J d 2 x d 2 y{^5 2 {x l - y[) - ^<$Vi - y'j] cos V^[<f>(x) + <j>{y) 
which leads to the flow equations for A and e. 



S' int , (5.13) 



11 



References 

[1] R.P.Feynman and J.A.Wheeler, Rev. Mod. Phys. 17 (1945) 157. 

[2] D.A.Eliezer and R.P.Woodard, Nucl. Phys. B 325 (1989) 389. 

[3] E.C.Marino, Phys. Lett. B 263 (1991) 63. Nucl. Phys. B 408 (1993) 551. 

[4] D.G.Barci, C.D.Fosco and L.Oxman, Phys. Lett. B 375 (1996) 267. D.G.Barci and 
L.Oxman, Mod. Phys. Lett. A 12 (1997) 493. 

[5] I. V. Krive et al, Phys. Rev. B 52 10865 (1995). 

[6] C.Naon, M.C.von.Reichenbach and M.L. Trobo, Nucl. Phys. B 435 [fs] 567 
(1995);ibid. B 485 [fs] 665 (1997). 

M.Mamas, C.Naon , M.L. Trobo, Nucl. Phys. B 525 [FS], 721 (1998). 

V. I. Fernandez and C. M. Naon, Theor. Math. Phys. (Russia), 118, 385 (1999), 

|hep-th/9810048| 

[7] S.Iijima, Nature 354, 56, (1991). 

T.W.Ebbesen, Physics Today 49 (6), 26, (1996). 

[8] H.Otani, T.Ogawa, Phys. Rev. B 54, 4540, (1996). 

[9] D. Poilblanc et al, Phys. Rev. B 56, R1645, (1997). 

[10] G. Fano et al, Phys. Rev. B 60, 15654, (1999). 
S.Capponi et al, ibid. 61, 13410, (2000). 

[11] R.Egger, A.O.Gogolin, Phys. Rev. Lett. 79, 5082, (1997). 

[12] B.Dumoulin et al, Phys. Rev. Lett. 76, 1360, (1996). 

[13] V.I.Fernandez, A.Iucci and C.M.Naon, Nucl. Phys.[FS]B 636, 514 (2002). 
|hep-th/0112057| 

[14] W.Thirring, Ann. of Phys. 3 91, (1958). 

B.Klaiber, Lectures in Theoretical Physics, Boulder, 1967, Vol 10A, Ed.A.Barut and 
W.Brittin, New York, Gordon and Breach, 141. 
V.E.Korepin, Teor. Mat. Fiz. (Russia) ,41 169, (1979). 

[15] J.Solyom, Adv. Phys. 28, 209 (1979). 

[16] A.Luther and V.J.Emery, Phys. Rev. Lett. 33, 589 (1974). 

R.Heidenreich, B.Schroer, R.Seiler and D.A.Uhlenbrock, Phys. Lett. A 54, 119 
(1975). 



12 



[17] V.J.Emery, A.Luther and I.Peschel, Phys. Rev. B 13, 1272, (1976). 



[18] Y. Saito, Z. Phys. B 32, 75, (1978). 

M.P.A. Fisher and W. Zwerger, Phys. Rev. B 32, 6190, (1985). 
M. Stevenson, Phys. Rev. D 32, 1389, (1985). 

R. Ingermanson, Nucl. Phys. B 266, 620 (1986). 

A.Iucci and C.M.Naon, J. Phys. A B 35, 8005 (2002). |hep^th/0209019| 

[19] R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 



H. J. Schulz, G. Cuniberti and P. Pieri, cond-mat/9807366 
S. Rao and D. Sen, |cond-mat/0005492| 

[20] S.Tomonaga, Prog. Theor. Phys. 5, 544 (1950). 
J.Luttinger, J. Math. Phys. 4, 1154 (1963). 
E.Lieb and D.Mattis, J. Math. Phys. 6, 304 (1965). 



[21] A. Iucci, K. Li and C.M.Naon, Nucl. Phys.[FS]B 601, 607 (2001)|h^th/0012269 

[22] Al. B. Zamolodchikov, Int. J. Mod. Phys. A 10, 1125, (1995). 

[23] See for instance, U. Ziilicke, Phys. Rev. Lett. 83, 5330, (1999); A. Rosch, 
cond-mat/0104181, to appear in Advances in Solid State Physics 41. 



13 



Figure captions 

Figure 1: Relative error Am/m as function of the stiffness constant, for the local case. 
The dashed line corresponds to the SCHA whereas the solid line shows the RG result. 

Figure 2: Phase diagram in the e — A plane, for K = .25. The dashed line is the 
critical line A = K e. 

Figure 3: Numerical solution of the gap equation for Ao = .5 and K = .25. The gap 
m decreases for increasing eo- 

Figure 4: Numerical solutions of the gap equations for A = .5, K = .25, and < 
e < 0.1. The dashed line corresponds to the SCHA whereas the solid line shows the RG 
result. The unit function is included to allow comparison with the local case. 
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